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We determine the solvent mediated contribution to the effective potentials for model colloidal or 
nano- particles dispersed in a binary solvent that exhibits fluid-fluid phase separation. The interac- 
tions between the solvent particles are taken to be purely repulsive point Yukawa pair potentials. 
Using a simple density functional theory we calculate the density profiles of both solvent species in 
the presence of the 'colloids', which are treated as external potentials, and determine the solvent 
mediated (SM) potentials. Specifically, we calculate SM potentials between (i) two colloids, (ii) 
a colloid and a planar fluid-fluid interface, and (iii) a colloid and a planar wall with an adsorbed 
wetting film. We consider three different types of colloidal particles: colloid A which prefers the 
bulk solvent phase rich in species 2, colloid C which prefers the solvent phase rich in species 1, and 
'neutral' colloid B which has no strong preference for either phase, i.e. the free energies to insert 
the colloid into either of the coexisting bulk phases are almost equal. When a colloid which has a 
preference for one of the two solvent phases is inserted into the disfavored phase at statepoints close 
to coexistence a thick adsorbed 'wetting' film of the preferred phase may form around the colloids. 
The presence of the adsorbed film has a profound influence on the form of the SM potentials. In 
case (i) reducing the separation between the two colloids of type A leads to a bridging transition 
whereby the two adsorbed films connect abruptly and form a single fluid bridge. The SM potential is 
strongly attractive in the bridged configuration. A similar phenomenon occurs in case (iii) whereby 
the thick adsorbed film on colloid A and that at the planar wall, which prefers the same phase as 
colloid A, connect as the separation between the colloid and the wall is reduced. In both cases the 
bridging transition is accompanied, in this mean-field treatment, by a discontinuity of the SM force. 
On the other hand, for the same wall, and a colloid of type C, the SM potential is strongly repulsive 
at small separations. For case (ii), inserting a single colloidal particle near the planar fluid-fluid 
interface of the solvent, the density profiles of the solvent show that the interface distortion depends 
strongly on the nature of the colloid-solvent interactions. When the interface disconnects from the 
colloid there is, once again, a discontinuity in the SM force. 

PACS numbers: 05.20.Jj, 82.70.Dd, 61.20.Gy, 68.08.Bc 



I. INTRODUCTION 

The effective forces between large bodies in a solvent, 
for example between suspended particles or between par- 
ticles and container walls, are composed of contributions 
from the 'direct' interactions, e.g. Coulomb and disper- 
sion forces and also from 'indirect' or solvent mediated 
interactions arising from the presence of the solvent. A 
significant part of colloid science consists of determin- 
ing, and often tailoring, the effective interactions between 
colloidal particlesP^l In this paper we focus on a par- 
ticular class of solvent induced interactions which arise 
from the adsorption of liquid films around large bodies. 
We consider large bodies suspended in a bulk solvent 
that exhibits coexisting fluid phases. The existence and 
thickness of the adsorbed 'wetting' films depend strongly 
on the statepoint of the solvent. By changing the sol- 
vent temperature or concentration such that the fluid ap- 
proaches coexistence, the thickness of the adsorbed film 
can increase and when two bodies with adsorbed films 
come sufficiently close these films can join together to 
form a fluid bridge. Such a mechanism generates strong 
attractions between the bodies that often lead to aggre- 
gation phenomena. Aggregation may be reversed by re- 



versing the temperature or concentration change - see 
Ref P for an overview of the purported role of bridging in 
colloidal flocculation. 

In the last decade or so interest has also grown in un- 
derstanding the effective interactions between colloidal 
and nano- particles and interfaces, and between pairs of 
particles adsorbed at an interface P Colloidal particles are 
often attracted to fluid interfaces and can st abilize emul- 
sions in a similar way to surfactants ISIS12I The adsorp- 
tion of colloidal and nano-scale particles at interfaces is 
important for a number of industrial processes includ- 
ing foams, lubrication, adhesion, and stabilizing cmul- 
s j ons JUsE The self assembly of nanoparticles and colloids 
adsorbed at fluid interfaces^' has been applied in many 
ways including creating solid 'capsule' structures.^ Par- 
ticles at interfaces are also used in fundamental studies 
of the phase behaviour and physical properties of (quasi) 
two-dimensional fluids and crystals)- 11 * 12 * 13 * 1 ^ Colloidal 
nanoparticles may also be used as building blocks for 
materials with specific mechanical, optical and magnetic 
properties.^ 

Here we investigate the effective interactions between 
'colloidal' particles and fluid interfaces using a micro- 
scopic treatment of a model system where the fluid in 
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which the 'colloids' are dispersed, henceforth referred 
to as the 'solvent', is composed of a binary mixture of 
soft-core particles which are much smaller than the 'col- 
loids'. These 'solvent' particles interact via purely repul- 
sive point Yukawa potentials which constitute a simple 
model f or the interactions between charged particles - 
see Refs.L!iy2G3 anc ] references therein. We showed previ- 
ously^ that for certain choices of mixing parameters this 
model solvent exhibits fluid-fluid phase separation and we 
focus our present study on state-points where the solvent 
is at, or close to, fluid-fluid coexistence. Our aim is to in- 
vestigate the effective solvent mediated (SM) interaction 
between large bodies ('colloids') in the presence of the 
adsorbed films. In the following we consider three phys- 
ical scenarios: (i) the effective interaction between two 
colloids of the same type immersed in the bulk solvent; 
(ii) that between a single colloid and the planar inter- 
face between two coexisting fluid phases; and (iii) that 
between a single colloid and a planar wall which is cov- 
ered by a thick adsorbed film of solvent. In each case we 
determine the equilibrium density profiles of the binary 
solvent in the presence of the colloids and from these pro- 
files we calculate the SM interactions. We use a simple 
microscopic density functional theory (DFT), shown to 
be reliable for bulk fluids, to calculate the density pro- 
files and grand potential of the solvent in the external 
potential of the fixed colloids. We describe the interac- 
tion between the colloids and the solvent particles by a 
hard-core potential plus a purely repulsive Yukawa tail 
that models crudely the effective interaction between a 
hard body carrying a surface charge and solvent particles 
with the same sign. Although our system may be con- 
sidered to be a toy model for charged colloidal particles, 
we believe that our results should provide insight into 
the behavior of any system of colloidal or nano- particles 
that are dispersed in a solvent exhibiting fluid-fluid phase 
separation. 

In the previous stud}^ we showed that the pair corre- 
lation functions of the bulk fluid at states close to coexis- 
tence are reasonably well described by a simple random 
phase approximation (RPA). Subsequently, using a den- 
sity functional theory that generates the RPA, we showed 
that the model fluid wets completely a hard-wall and 
exhibits a pre-wetting transition slightly away from co- 
existence.^ Furthermore, by adding a repulsive tail to 
the wall potential we found that the location and ex- 
tent of the pre-wetting line in the phase diagram can be 
changed. Although the model solvent is very simple, its 
bulk and interfacial phase behavior mimics that in more 
sophisticated models and has the important advantage 
that it can be studied using the computationally inex- 
pensive RPA based DFT. 

We consider three different colloids that vary in their 
affinity for one of the coexisting bulk fluid phases. Col- 
loid type A has a very strong preference for the phase 
rich in species 2, colloid B is almost 'neutral' but has a 
very weak preference for the phase rich in species 2, and 
colloid C has a strong preference for the phase rich in 



species 1. When colloids A and C are inserted into their 
unfavored phase they adsorb a thick layer of their pre- 
ferred phase. We find that the propensity of an isolated 
colloid immersed in the bulk phase to adsorb thick films 
dominates the behavior of the solvent when the colloids 
are brought together, or brought to interfaces, and there- 
fore determines to a large extent the nature of the SM 
interactions between bodies. In particular we find that 
for two colloids in bulk the presence of thick adsorbed 
films leads to a long-ranged and strongly attractive SM 
interaction: as the two colloids are brought together the 
adsorbed films join together abruptly to form a bridge 
between the colloids thereby minimizing the interfacial 
contribution to the grand potential. If in isolation the 
colloids do not adsorb thick films then the SM poten- 
tial is much shorter ranged but is still attractive due to 
depletion effects. For a colloid that strongly prefers to 
be in one of the phases the effective interaction with the 
fluid-fluid interface is highly asymmetric with the global 
minimum of the grand potential, that determines the col- 
loid's equilibrium position, occurring in bulk, far away 
from the interface. For the 'neutral' colloid the global 
minimum of the SM interaction is in the center of the 
interface where the colloid intersects the interface. For 
the third scenario, a single colloid near a planar wall, we 
find that if the colloid prefers the phase adsorbed at the 
wall then there is a long-ranged, attractive SM interac- 
tion between the wall and the colloid. If the colloid does 
not strongly prefer either phase it seeks a location that 
intersects the interface between the thick adsorbed film 
and the bulk fluid and the interaction is attractive but 
not as strong, nor as long-ranged, as the previous case. 
On the other hand, if the colloid prefers the bulk phase 
then the SM potential will be repulsive, and relatively 
short-ranged. 

The use of DFT for studying the wetting behavior of a 
solvent in the presence of an isolated big spherical parti- 
cle that exerts a (spherical) external potential on the sol- 
vent atoms is, of course, well established see e.g! 18 * 19 * ^ 
and references therein. However, the behavior of the sol- 
vent in the presence of two fixed particles or a particle 
and a wall is not as widely studied. These more complex 
geometries generally require considerably more computa- 
tional effort. In the following we mention several different 
approaches that have been used to investigate these prob- 
lems using DFT. One of earliest studies of the SM inter- 
action between two spherical particles exerting dispersion 
forces made use of an interface potential approach where 
the solvent density profile between liquid a nd ga s phases 
is modeled by a sharp-kink approximationPi!221 gy so i v _ 
ing for the position of the interface around a pair of col- 
loids at decreasing separations an abrupt crossover from 
unconnected to connected liq uid fil ms, termed capillary 
condensation by the authors o f 21 * 22 !, was found. In RefP^ 
a more sophisticated free-energy functional was employed 
but the sharp-kink interface approximation was retained. 
The authors investigated the morphological transition in 
some detail and also compared the SM interaction poten- 
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tial to the bare interaction arising from dispersion forces 
between the pair of spheres. 

An alternative DFT approach developed by Roth et 
aP^l makes use of the potential distribution theorem 
which provides a formally exact expression for the SM 
potential in terms of the one-body direct correlation func- 
tion between the inserted particle and the solvent. This 
in turn depends only on the inhomogencous solvent den- 
sity profile around a single, isolated (spherical) particle 
before the second particle is inserted. Given a suitable 
density functional, the direct correlation function and 
the spherically symmetric density profile of the solvent 
around a single colloid are easily determined and it is 
straight forward to calculate the SM potential between 
two large particles. The advantage of this approach is 
that one must only calculate a solvent density profile 
that has spherical symmetry. The disadvantage is that 
one requires a DFT reliable for highly asymmetric mix- 
tures. This so called insertion method has been used to 
calculate the SM potenti al for v arious types of hard and 
soft-core model particles! 24 * 25 * 26 * 

In Ref!^3 Archer et al. compared the results of the in- 
sertion method to those of a 'brute-force' DFT approach 
for calculating the SM potential between two large Gaus- 
sian soft-core particles immersed in a bulk binary sol- 
vent of smaller soft-core particles. Note that we refer to 
the procedure that calculates the solvent density profiles 
around two fixed colloids treated as external potentials 
as the 'brute-force' method, in contrast to the insertion 
methocPS that requires only a calculation of the solvent 
density profiles around a single fixed colloid. The au- 
thors found that for certain solvent statepoints close to 
fluid-fluid coexistence the large particles could adsorb a 
thick film of the wetting 'phase'. Using the brute force 
approach, it was shown that on reducing the separation 
between the two large particles the adsorbed films would 
connect abruptly forming a fluid bridge. Such a bridging 
transition leads, at mean-field level, to a discontinuity 
in the SM force between two large particles and in the 
bridged phase the potential is long-ranged and strongly 
attractive. Results from the insertion method, using the 
same RPA based DFT but now for a ternary mixture, 
shown to be accurate for solvent statepoints away from 
the binodal, could not account for the bridging transi- 
tion. However, the insertion method does predict, albeit 
qualitatively, the occurrence of long-ranged SM poten- 
tials when the thick adsorbed films are present.^ 

A number of other studies have used brute-force DFT 
methods to calculate the SM potential. Stark et al. 28 
used Landau-de Gennes theory to investigate capillary 
bridging in the case of two large hard spherical colloids 
immersed in a bulk isotropic liquid crystal host at state- 
points close to the isotropic-nematic phase boundary. 
GrodorP^ studied a model colloid-polymer mixture ad- 
sorbed between two large colloidal spheres. Andricnko 
et alW^ used Landau theory to calculate density profiles 
and force distance curves for a solvent adsorbed between 
a large sphere and a planar wall. They found a first- 



order capillary bridging transition and investigated its 
depe ndence on the radius of the sphere. Cheung and 
Aii er (3!]32i33| nave usec j an Onsager (second virial) DFT 
for hard rods to investigate the effective interactions be- 
tween two cylindrical colloids, and between a cylindrical 
colloid and a wall, in a liquid-crystal host. At suitable 
statepoints in the bulk isotropic phase, regions of nematic 
order can form at the wall and around the cylinders. The 
adsorbed nematic regions can then form a bridge, in an 
analogous fashion to the above examples, giving rise to 
strongly attractive interactions. 

We are not aware of studies that use DFT to study the 
interaction between a large particle and a fluid interface 
although we are aware of studies that use self-consistent 
field-theory (SCFT) to investigate the interactions be- 
tween particles and interfaces in self-organized block co- 
polymer structures! 34 * 35 ^ 

There is a growing body of computer simulation in- 
vestigations concerned with the adsorption of individ- 
ual s pherical, and non-spherical, particles at fluid inter- 
faces pEM3MMEo] j n general the size of the nanoparticles 
is of the same order as the fluid particles and thick ad- 
sorbed layers are not present or are not considered. Much 
is made of li ne ten sion contributions to the free energy 
of immersion! 36 * 40 * 

Our paper proceeds as follows: In Sec. |TT] we recall 
some of the basic background to the thermodynamics 
of solvation, colloid-colloid and colloid-interface interac- 
tions. In Sec. |III| we outline the model system and in Sec. 
|IV| we describe our density functional theory approach. 
Sec. [V] describes the results of this investigation. Finally 
in Sec. I VII we discuss the results and make some con- 
cluding remarks. Note that we frequently use the term 
'colloid' for the particle or particles inserted into the sol- 
vent. However, the inserted particle has an effective ra- 
dius of 4.5A -1 whereas that of the solvent particles is 
about 0.5A -1 . Such length scales are much more appro- 
priate to those associated with inserting nanoparticles 
into solvents. 



II. BACKGROUND: SOME KEY CONCEPTS 

A. Thermodynamics of solvation 

Before describing our microscopic DFT approach for 
calculating SM potentials, we first recall some of the ther- 
modynamics relevant to a colloid (spherical particle) im- 
mersed in a fluid. We consider a grand canonical system 
enclosed within a volume V containing several different 
species of solvent particles, that is coupled to a reservoir 
which fixes the temperature T and the set of chemical 
potentials {[ii} for the different solvent species i in the 
system. Later in this paper we consider the specific case 
of a solvent composed of two different species of particles 
(i.e. i = 1 or 2), but for now we leave the number of com- 
ponents undefined. At equilibrium, the grand potential 
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of the uniform solvent is 

il(V,T,{iH}) = -PV, (1) 

where P is the pressure and we have assumed that there 
are no interfacial contributions. If we insert a spherical 
particle (colloid) with radius R, the grand potential is 

n(V,T,{m},R) = -P(V- ^nR 3 ) + 4nR 2 lpa (R), (2) 

where (V - ±irR 3 ) is now the volume that is accessible 
to the solvent (which we denote phase a) and "f pa (R) 
is the surface tension (surface excess grand potential per 
unit area) between the colloidal particle p and phase 
Note that when R is large, and in the absence of wetting 
films: 

lpa(R) = l pa (oo) (l - ^ H ^ , (3) 

where 6 is the Tolman length, i.e. a microsco pic lengt h 
scale of order the radius of the solvent particles! 19 * 20 * 41 ! If 
we subtract Eq. from Eq. ^ , we obtain the following 
expression for the excess grand potential for inserting the 
particle into the system (phase a): 

^a x = \*R 3 P + AixR 2 lpa {R). (4) 

We now consider the case when the chemical potentials 
{/Ltj} and the temperature T are such that we observe two 
phase coexistence between phase a and a second phase 
which we denote (3. In this case, the grand potential of 
the system is 

tt(y,T,{i H }) = -PV + j a pA, (5) 

where 7 Q/ 3 is the surface tension for the (planar) interface 
between phases a and (3 and A is the area of this interface. 
If we insert the colloid into phase a, far away from the 
interface, then the excess grand potential is still given by 
Eq. Q, but if we insert the colloid into phase /3, again 
at a point far away from the interface, then the excess 
grand potential is now 

flf = ^irR 3 P + 4nR 2 lpf3 (R), (6) 

where r y p p(R) is the surface tension between the colloidal 
particle p and phase (3. 

When the two surface tensions 7 pct and 7 P/ 3 are equal, 
then the colloid does not favor one phase over the other. 
However, this is a special case and more generally the 
surface of the colloid will favor one phase over the other. 
Suppose it favors phase a, i.e. 7 pQ < 7 P( g. There is 
an upper bound on the value 7 P( 3 can take, since in the 
case when the surface of the colloid strongly favors the 
a phase, on inserting the colloid into the coexisting (3 
phase, a film of the favored a phase forms around the 
colloid. Thus, the upper bound on 7 Pi a is given by 7 P) 9 < 



"/ pa +"fai3- The equality occurs in the limit that the radius 
R — > oo, and when the surface is completely wet by phase 
a. For a colloid of radius R having this kind of surface, 
we find that the excess grand potential for inserting it 
into the f3 phase is: 

nf^n e a x + 4Tr(R + i) 2 la p(R + i), (7) 

where f2® x is given by Eq. Q , Z is the thickness of the film 
of phase a covering the surface of the colloid and 7 Q| a(i?) 
is the surface tension of the spherical fluid-fluid interface 
between the wetting a phase and the bulk (3 phase. Note 
that in Eq. we have omitted the interaction between 
the two interfaces pa and a(3. This result also holds when 
the solvent (3 phase is not exactly at coexistence: chang- 
ing the chemical potentials so that the mixture is (not 
too far) off coexistence simply reduces the wetting film 
thickness Z. However, there is an additional contribution 
proportional to Z a rising from the fact that the a phase 
is metastablePMl 



B. Solvent mediated potentials in a bulk fluid 

We now consider inserting two colloids into the bulk 
fluid. When the pair of colloids are far apart (separated 
by a distance h — > oo), then the insertion free energy 
fl cx (h — > oo) is simply twice the result in either Eq. Q or 
Eq. (|6j, depending into which phase a or (3 the colloids 
are inserted. However, in the case when we insert the 
two colloids separated by a finite distance, the insertion 
free energy f2 ex (/i) now depends on the distance between 
the centres of the colloids h and also on the bulk phase 
into which the colloids are inserted. The SM potential is 
defined as follows: 

W(h) =n cx (h)-n ex (h^oo), (8) 

i.e. it is the difference between the grand potential when 
the colloids are far from one another (h — > oo) and when 
they are separated a distance h. The total effective po- 
tential between the two colloids is U(h) + W(h) where 
U(h) is the 'bare' or direct interaction. 

When the colloids are inserted into the preferred a 
phase W(h) is fairly short ranged. Due to depletion 
forces^ W{h) is typically attractive when the colloids are 
close to contact. Generally W(h) has a range determined 
by the bulk correlation length in the solvent. However, 
if we insert the two colloids into phase (3 close to or at 
bulk coexistence, where both colloids are surrounded by 
a thick adsorbed wetting film of phase a, then W(h) can 
become much greater in magnitude, much longer ranged 
and strongly dependent on the conformations of the thick 
films surrounding the colloids. When the colloids are far 
apart the adsorbed films do not interact, i.e. when h is 
large W(h) ~ 0. As the colloids are brought closer to- 
gether, the surrounding wetting films begin to interact 
and if the colloids are sufficiently close the system can 
minimize the interfacial area between the wetting a and 
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bulk (3 phases (and therefore minimize the grand poten- 
tial) by bridging the gap between the colloids and creat- 
ing a single interface. Without explicitly solving for the 
minimal interfacial shape as a function of h (see Eq. ( |25[ ) 
below for an approximate analytic expression for S(h)), 
we may write the excess grand potential for inserting the 
two colloids into the (i phase as follows [c.f. Eq. Q]: 

Qf(h) ~ 20™ + Sn(R + l) 2 -/ a p(R + l)S(h), (9) 

where S(h) is the ratio of total surface area of the ad- 
sorbed wetting film(s) when the particles are a distance 
h apart to that when they are infinitely far apart. In Eq. 
([9]) we have assumed that the average surface tension (ex- 
cess grand potential per unit area) of the a-(3 interface is 
equal to r y a /3(R + I) for all surface conformations. When 
the adsorbed films do not interact, S(h) = 1. However, 
when the colloids become close enough for the films to 
connect and form a bridge, then 0.5 < S(h) < 1. Bridge 
formation leads to a discontinuity in the first derivative 
of S(h) at h = h br . For h < h br , S(h) decreases smoothly 
as the colloids are brought to contact. Using Eq. ([9]), we 
can approximate the SM potential as follows: 



W{h) ~ 87r(R + l) 2 laf j(R + l)(S(h) - 1) 
~ -8nR 2 lafJ [l - S(h)\, 



(10) 



where we have assumed in the second line that "/ a /3(R + 
I) ~ 7 Q/ 3(oo) = "/ a /3, the surface tension of the planar a- 
f3 interface, and we have used the fact that R^> I. This 
means that it is the quantity A = R 2 "/ a /3 which deter- 
mines the value of W{h) when the colloids are surrounded 
by thick wetting films and close to contact, h ~ h c . Since 
typically the interfacial tension j a p ~ ksT/a 2 , where fee 
is Boltzmann's constant and a is a molecular length scale 
(the size of the solvent atoms/molecules), we find that 
A ~ fcsT(i?/cr) 2 . Thus for even modest sized colloids of 
radius lOnm, since a ~0.2nm for a typical solvent, the 
gain in free energy in bringing the colloids close to contact 
can be large; in this case, W(h c ) ~ — 10 3 to — 10 4 fcsT. 
It is for this reason that one should expect colloids to 
aggregate strongly when the solvent is near to phase co- 
existence in cases where the surface of the colloids is such 
that they are wet by the coexisting phase. 



C. Interaction between a colloid and an interface 

For colloids dispersed in a solvent which is at bulk 
phase coexistence (we continue to denote the two coex- 
isting phases a and p) one can consider the effective in- 
teraction between the colloid and the interface between 
the two phases. One generally finds that large colloidal 
particles are strongly attracted to the interface and be- 
come strongly bound to it. The reason for this is that 
if one places the colloid within the interface, so that the 
interface intersects the colloid, then the area of the in- 
terface between fluid phases a and (3 is reduced by an 



amount ~ irR 2 and so there is a change in the free en- 
ergy 



An ~ -7ri? 2 7 Q/3 



(11) 



on moving the colloid from the bulk fluid into the in- 
terface. Again we find it is the quantity A — R 2 "/ a /3 
which determines the magnitude of Af2. For colloids of 
radius lOnm, the change in grand potential Af2 ~ — 10 3 
to -10 4 /c B T. 

Going beyond this very rough approximation, we can 
consider a single large spherical colloid embedded within 
a non-deformable planar fluid-fluid interface: the excess 
grand potential as a function of the variable z — z Q /R, 
where zo is the distance of the centr e of t he colloid from 
the plane of the interface, is given by^ :J 



2nR l 



(1 - Z )(lp/3 ~ Ipa) 



■i(i -z 2 )j af3 + ~yj 



(12) 



when the colloid is within the interface, —1<zq<1. 
f2® x is given by Eq. Q and r is the line tension. The line 
tension was originally introduced by Gibbs to describe 
the excess free energy associated with the line where three 
phases meet P»Q This quantity may be positive or negative 
and experimental measurements of the line tension have 
proved to be difficult. 43 For large R the contribution from 
the line tension to Eq. (12) will be small compared to 
the contributions from the surface tensions. However, 
for nanoparticles line tension effects can be significant; 
see RefP, and references therein. 

By considering the interface to be a deformable mem- 
brane, one can go beyond Eq. ( [12] ) and include the con- 
tribution s to O ex from the bending and stretching of the 
interface This approach has also been extended to de- 
termine the effective potential between two colloids that 
are within an interface However, this interesting prob- 
lem is beyond the scope of the present study, though we 
return briefly to it in Sec. |VI| 



III. MODEL SYSTEM 



Solvent 



The solvent is a binary mixture of particles that inter- 
act via purely repulsive, point Yukawa pair potentials 



eMij exp(— Ar) 



4tt 



Ar 



(13) 



where the parameter e > sets the energy scale, My is 
the interaction magnitude between species i and j, and 
A is an inverse length scale. In RefP^ we showed that 
when M\2 > \[M^[M^2 and the bulk fluid total number 
density p b = pi + p 2 is sufficiently high, then the system 
phase separates into two fluid phases, one phase rich in 
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FIG. 1: Bulk phase diagram of the binary point Yukawa fluid 
(our solvent) with parameters Mil = 1, M\2 = 2.2, M22 = 4 
and temperature ksT/e = 1. p b is the total bulk density 
and X2 is the concentration of species 2. The de-mixed re- 
gion is bounded by the binodal (solid line) and meets the 
spinodal (dashed line) at the critical point, indicated by o. 
The straight (dotted) tie-lines connect coexisting state-points 
with pressures f3\~ 3 P = 150 to 500, in increments of 50, and 
then from 500 to 1900 in increments of 100 (from bottom to 
top). The crosses mark three particular state-points: the co- 
existing pair X and Y at the pressure /3A~ 3 P = 233, which 
are joined by the dash-dot tie-line and the point W that lies 
very close to coexistence at p b \~ 3 — 24.0. Further details 
concerning these three state points, at which we display var- 
ious results below, are given in Table [I] The inset displays a 
magnification of the region around points W and X. 



particles of species 1 and the other phase rich in par- 
ticles of species 2. The phase diagram for the mixture 
with parameters M\\ = 1, M22 — 4, M12 = 2.2 and tem- 
perature kgT/e — 1, calculated using the random phase 
approximation (RP A) for the Helmholtz free energy (see 
RefPl and Sec. 
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below for further details) is displayed 
in Fig. [T] Clearly this system exhibits fluid-fluid phase 
separation. Indeed this is the main reason for its selec- 
tion as model solvent in this study. Another important 
reason for studying this model fluid is that in Ref we 
showed that when this system is at coexistence there can 
be complete wetting of the planar hard-wall by the phase 
rich in species 2 accompanied out of coexistence by a first- 
order surface (pre-wetting) phase transition from a thin 
to a thick adsorbed wetting film. We showed that this 
was also the case for a planar hard wall potential aug- 
mented with repulsive exponential or Yukawa potentials. 
We found that the location and extent of the pre-wetting 
line depends strongly on the parameters of the wall po- 
tential.^ These properties allow us to use this micro- 
scopic model to investigate the influence of wetting films 
on solvent mediated potentials between colloidal parti- 
cles. 



B. Colloids 

We model the interaction of the colloids with the 
solvent particles i —1,2 via a hard-core plus repulsive 
Yukawa-tail pair potential: 



#°(r) 



, exp(-A(r-.R)) 
v fcu * A(r-fl) 



< r < R 
R<r, 



(14) 



where R is the radius of the colloids and the parameter 
bi > determines the magnitude of the repulsive tail in- 
teraction with solvent particles of species i. The value of 
the ratio 61/62 is important for determining whether the 
surface of the colloid favors solvent particles of species 1 
or particles of species 2. In this study we choose both 61 
and 62 to be fairly large in order to 'soften' the bound- 
ary of the hard-core. If the potentials in Eq. ( 14 1 change 



rapidly over a short distance (which is the case when 
both 61 and 62 are small), and the potentials are de- 
fined on a Cartesian grid this can lead to numerical er- 
rors in the solvent density profiles, which are defined on 
the same Cartesian grid. All the results presented in this 
paper are for the case when the radius of the colloids 
R = 4A _1 . The effective radius of the solvent particles 
is roughly between 0.2A -1 - 0.5A -1 for the densities of 
interest here. We estimate this radius from the size of 
the 'correlation hole' in the radial distribution functions 
9ij( r ) ~ see RefP^for examples of gij{r) for this system. 
A radius R = 4A _1 is sufficiently large for the colloid to 
be covered by a thick adsorbed 'wetting' film when the 
ratio 61/62 is such that the colloid favours one solvent 
species over the other and the solvent is at a state point 
near to coexistence. In this study, we set 61 = 10 and 
62 = 10, 20 and 30, thus defining three different types 
of colloids, A, B and C. We find that colloids of type A, 
with 62 = 10, have a strong preference for solvent par- 
ticles of the 'bigger' species 2. Colloids of type B, with 
62 = 20, have a very weak preference for species 2; B 
colloids are essentially 'neutral' colloids. Colloids of type 
C, with 62 = 30, have a strong preference for species 1, 
the 'smaller' solvent particles. Note that we do not define 
the 'bare' interaction potential U(h) between the colloids. 
This quantity does not enter into our calculation of the 
solvent mediated potential W(h). 



C. Wall-Solvent Interfaces 

We model the potentials exerted by a planar wall on 
the solvent particles as follows: 

I Oieexp(— Xz)/[Xz) z > 0, 

where the parameter aj > determines the magnitude 
of the repulsive tail interaction with solvent particles of 
species i. In a manner entirely analogous to the way the 
colloid potentials in Eq. (14 1 influence the solvent, the 
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value of ratio a\ / ai is important for determining whether 
the wall favors solvent particles of species 1 or solvent 
particles of species 2, i.e. the precise values of the pa- 
rameters a\ and a-i determine whether there is a wetting 
film of the coexisting phase adsorbed at the wall when 
the bulk phase is brought to coexistence and also deter- 
mine the location and extent in the phase diagram of the 
pre-wetting phase transition lineP^ In this study we set 
ax = 02 = 1. For these values we find that there is a 
weak preference for particles of species 2 at the wall, so 
that when the bulk phase rich in particles of species 1 
is brought near to coexistence, we find that there is a 
thick wetting film of the phase rich in species 2 adsorbed 
at the wall. Any wetting transition occurs at densities 
above those shown in Fig. [T] These values of a! and 
ai are sufficiently large that both solvent fluid density 
profiles vary sufficiently slowly for there not to be any 
numerical errors due to the discretisation of the density 
profiles. 



IV. DENSITY FUNCTIONAL THEORY AND 
ITS IMPLEMENTATION 

Here we present a brief description of the DFT ap- 
proach that we use to calculate the solvent density pro- 
files around the colloid(s) and to calculate thermody- 
namic quantities such as f 2 ex . For a more complete ac- 
count of DFT sec e.g. For a fluid composed 
of v different species of particles one can construct a 
functional f2[{/3i}] of the set of one-body density profiles 
{pi(r)}, i — 1 . . . v. The minimum value of this func- 
tional is equal to the thermodynamic grand potential of 
the system, i.e. f2 = min(f2 [{/?$}]), and the set of density 
profiles which minimize the functional f2[{pi}] are the 
equilibrium fluid density profiles satisfying the following 
set of Eulcr-Lagrange equations: 



6to[{pi}] 



0, 



(16) 



which may be solved to obtain the set of equilibrium fluid 
density profiles. The grand potential functional may be 
written as: 



"[{ft}] = F[{ Pl }} dr Pi (r)( Mi 



Vi(r)), (17) 



where ^[{pi}] is the intrinsic Helmholtz free energy func- 
tional and Vi(r) is the external potential acting on par- 
ticles of species i. The intrinsic Helmholtz free energy 
functional may be separated into a sum of two contri- 
butions: ^[{p*}] = ^id[{Pi}] + F cx [{pi}} where the first 
term is the Helmholtz free energy of an ideal-gas 

F id [{ Pi }] =Yk B T [ dr ft (r)(ln(Af^(r)) - 1), (18) 



where Aj is the thermal de Broglie wavelength of particles 
of species i, and the second term is the excess contribu- 
tion due to the particle interactions. 

The solvent that we consider is a binary mixture 
[y = 2) of soft-core particles interacting via the pair po- 
tentials in Eq. (13 1. For this system, at the state points 
of interest, it is sufficient to employ a simple mean-field 
approximation for the excess part of the Helmholtz free 
energjff^l 

2 

FAiPr}] = \J2 J dr J drV*(r)p,(r')<M|r-r'|). 



»ii=i 



(19) 

This functional generates the RPA for the pair direct 



-1 12ll9l 



S 2 F^{{ Pl }] 
<5p;(r)<5pj(r') 



correlation functions: cy(|r — r' 

-^Qr-r'l), where (3 = {k B T) 

In order to simplify the computation of the density 
profiles it is convenient to reduce the range of the solvent 
pair potentials 4>ij{r) by cutting and shifting these, i.e. 
we replace the solvent pair potentials 4>ij(r) by </>y s (r), 
which arc defined as follows: 



C 8 W 



[<t>ij(r) - <f>ij(r c )] r<r c 







r > Tr 



(20) 



where 9 ~ 1 is a dimensionless factor. If we set 9 = 1, 
then the system with the cut and shifted (CAS) potential 
does not have the same bulk phase diagram as the fluid 
with the full range potential <j>ij(r) (one is shifted slightly 
with respect to the other). However, if we choose the 
value of 9 so that 



&rr 2 (t>ij(r) 



drr 2 C as (r), 



(21) 



then within the RPA, the bulk phase diagram for the fluid 
with the CAS potential is identical to the phase diagram 
of the system with the full range potential. This makes 
comparisons with existing results more straightforward. 
In the present study we set the cut-off distance r c — 
4A _1 . This is large enough so as to not alter significantly 
the correlation functions of the fluid. 

The colloids are treated as fixed external potentials, 
which must be added to any other external potentials, to 
give 



^(r) = ^ s (r) + ]TC(|r-r fe |) 

k 



(22) 



where V^ ys (r) is given by Eq. (15) (when no wall is 
present, V^ sys (r) = 0), </>f c (r) is given by Eq. (14 1 and 



is the location of the centre of colloid k. In the present 
study we only consider situations that have cylindrical 
symmetry, i.e. one colloid plus a planar interface or an 
isolated pair of colloids, and so it is natural to character- 
ize points in space by the coordinates (z,r), where z is 
the axial distance along the z-axis, which passes through 
the center of the colloid (s) and is perpendicular to the 
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plane of the interfaces (when interfaces are present), and 
r is the radial distance from the z-axis. We denote the 
location of the center of colloid k along the z-axis by Zk- 

The solvent density profiles are calculated numerically 
by solving Eqs. (16) self consistently using a simple Pi- 
card iteration algorithm. We use fast Fourier transforms 
to evaluate the convolution integrals in the axial direc- 
tion, but the radial convolutions are integrated directly. 
This allows us to use a grid with a smoothly varying dis- 
tance between grid-points in the radial direction. Close 
to the colloids we use a small grid spacing ~ 0.1A -1 in 
order to accurately determine the rapidly varying density 
profiles at these points, but further out, where the pro- 
files vary slowly, we use a larger grid spacing ~ 0.3A -1 in 
order to increase the computation box size without any 
loss of performance. 

In calculating the solvent density profiles around the 
colloid in the fluid-fluid interface it is necessary to fix the 
density profiles on the boundary of the system far away 
from the colloid along the r-axis to be equal to those in 
the unperturbed state, when the colloid is absent, with 
the interface fixed at a certain value of z: 



Pi(z,r = r L ) = pi(z,r -> oo), 



(23) 



where rj, is the size of the computation region. If this 
is not done, the entire interface simply translates to the 
value of z where the free energy is the global minimum 
rather than giving the constrained free energy minimum 
where the colloid is constrained to be a certain distance 
from the plane of the interface. For the particle in the 
fluid interface calculations we set ~ 70A -1 . For the 
case of a colloid at a wall, provided that the size of the 
computation region is sufficiently large, the density pro- 
files satisfy Eq. (23 1 automatically, since the location of 



the interface between the wetting film and the bulk fluid 
is determined by the properties of the wall and the chem- 
ical potentials. 

In order to calculate the surface tension (excess grand 
potential) 7 associated with the individual interfaces, and 
compare with results such as Eq. (10 1, we use the equi- 



librium fluid density profiles together with the following 
expression: 



1 = 



n + pv 



dr( W [{ ft (r)}] + P) 



(24) 



where A is the area of the interface, 
grand free energy density [see Eq. 
— aj[{pi(oo)}] is the bulk fluid pressure 



uMpi 



(r)}] is the 
and P = 



V. RESULTS OF CALCULATIONS 

A. Interfaces in the absence of colloids 

Before examining the solvent mediated interactions be- 
tween the colloid and the interfaces, we first describe 



System 


Labels 


p b \- 3 


X2 


/3 7 A" 2 


Thick film at wall 


W 


24.0 


0.0112 


-116 




X 


20.0 


0.0267 




Fluid-fluid interface 


5.46 




Y 


11.8 


0.766 





TABLE I: Two systems are considered: (i) state point W 
in Fig. [TJ where the concentration X2 is slightly less than the 
value at coexistence a^.cocx = 0.01132 and the planar wall is 
covered by a thick adsorbed wetting film rich in particles of 
species 2, and (ii) the fluid-fluid interface between co-existing 
phases at state points X and Y in Fig.[T] 7 is the excess grand 
potential per unit area (surface tension) for these interfaces 
and is calculated using Eq. (24 1. 



briefly the properties of the solvent at the interfaces with- 
out any colloids present. Full details of the solvent inter- 
facial behavior can be found in RefP^ here we merely 
recall the properties relevant to the present study. 

In Fig. we display the solvent density profiles at 
the planar fluid-fluid interface between the coexisting 
state points X and Y, which have a bulk fluid pressure 
/3PA~ 3 = 233. The values of the total density and con- 
centration at each of these two state-points are listed in 
Table [T] The density profiles are monotonic functions of 
z and the (10%-90%) width of the interface is 4.2A" 1 . 
These density profiles are typical of a binary system that 
is at a state point not too far removed from the bulk fluid 
critical point. 

For the bulk fluid at state point W (see Fig. [T]) in 
contact with a planar wall interacting with the solvent 
via the potentials in Eq. (151, with eii = (12 = F we 



find that the wall is covered by a thick adsorbed wetting 
film rich in particles of species 2 - see the density profiles 
displayed in Fig. [2jb) . The thickness of the wetting film is 
determined by the separation of the bulk state point from 
the binodal; the closer the state point is to the binodal, 
the thicker the wetting film.EO State point W is chosen 
because here the wall is covered by a relatively thick film 
with a well-defined interface boundary; recall that as one 
approaches the critical point, the width of the interface 
between the adsorbed film and the bulk fluid grows and 
ultimately diverges, as does the thickness of the adsorbed 
film. In Fig. [2]d we also see that there is a marked peak 
in the density profile for species 2, corresponding to a 
strongly adsorbed layer of particles at the wall. Beyond 
this layer, the thick wetting film extends out a distance 
~ 5A _1 from the wall. 



B. A single colloid in the bulk solvent 

We now describe the properties of the solvent around 
a single colloid immersed in a bulk phase that is at coex- 
istence. Using the solvent-colloid pair potentials and the 
parameters defined in Sec. |IIIB , we calculate the density 
profiles of the solvent around a single colloid at the co- 
existing state-points X and Y. Our results are displayed 
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30 



Statepoint X 



Statepoint Y 



FIG. 2: Solvent density profiles pi{z) for two inhomogeneous 
systems: a) The planar fluid-fluid interface between coexist- 
ing phases at state-points X (rich in species 1) and Y (rich 
in species 2) in Fig. [I] b) The solvent density profiles at a 
planar wall with potentials given by Eq. (15 1, for the bulk 
fluid at state point W in Fig. [I] These profiles show that 
there is a thick wetting film rich in solvent particles of species 
2 adsorbed at the wall. Note that the density profiles of both 
species vanish at the wall. 



in Fig. [3] At both state points, we see in Fig. [3ja) and 
(b) that the surface of colloid A strongly prefers solvent 
particles of species 2; the density of species 1 particles 
at the surface of the colloid is much lower than that of 
species 2. This leads to the adsorption of a thick film of 
the coexisting phase rich in species 2 when colloid A is 
inserted in the bulk fluid at state-point X. Recall that the 
thickness l(R) of this wetting film is determined by the 
curvature (i.e. radius R) of the surface of the colloid 
In the limit that the colloid radius R — > 00, the thick- 
ness of the wetting film of the coexisting phase diverges 
(becomes macroscopically large): for short-ranged forces 
of the type considered here I ~ ln(i?). 20 Colloid B has a 
weak preference for particles of species 2, but in Fig.|3jc) 
and (d) we see only a small increase in the density of both 
species near the surface of colloid B; there are no thick 
wetting films. We see in Fig.[3je) and (f) that the behav- 
ior of the solvent around colloid C is essentially the oppo- 
site of the behavior around colloid A: particles of species 1 
are strongly adsorbed at the surface, so that when colloid 
C is immersed at state-point Y, there is a thick adsorbed 
wetting film of the coexisting phase rich in species 1 par- 
ticles around the colloid. From the density profiles we 
estimate an effective colloid radius, R* , by taking the dis- 
tance where the density of the strongly adsorbed species 
is 10% of its maximum value. For all colloids, in both 
phases, this occurs at approximately r = R* = 4.5A -1 . 
We also compute the grand potential ft for colloids A, B 
and C within the bulk solvent phases at points X and Y. 
This allows us to quantify the preference of the colloids 
to be in either bulk phase by calculating the difference 
in grand potentials, Ail — il Y — il x . Colloid A has a 
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FIG. 3: The solvent density profiles pi(r) around colloids A, 
B and C of radius R = 4A _1 at the state-points X and Y 
in Fig. [T] (see also Table [l| . r is the distance from the centre 
of the colloid. See the text for a discussion of these density 
profiles. 



very large negative value, Ail = —4870kgT, indicating 
that it strongly prefers to be in bulk phase Y. Colloid B 
has small negative value, Ail = — 403/csT, indicating a 
slight preference for phase Y. The grand potential differ- 
ence for colloid C is Ail = 4610kgT. This large positive 
value indicates that colloid C had a strong preference for 
phase X. Note that by phase X (Y) we mean the bulk 
phase at state-point X (Y). 



C. A pair of colloids in the bulk solvent 

We now describe the behavior of the solvent for two 
colloids separated by a distance h in the bulk fluid at 
statepoint X. The SM potential between them, W(h), 
is denned in Eq. In Fig. [4] we display the solvent 

density profiles around two colloids of type A at vari- 
ous separations h between the centers of the colloids and 
in Fig. [6] we display the SM potential W(h) obtained 
from these density profiles. As discussed in the previ- 
ous subsection, we find that the colloids are covered by 
a thick adsorbed wetting film that is rich in particles of 
species 2. For large separations h > 19A -1 the films do 
not interact with each other, see Fig. [4] and the potential 
W(h) — 0. As the separation h is decreased, the wet- 
ting films around the two colloids begin to influence one 
another and when the colloids are brought to a separa- 
tion h = /ibr = 18.7A -1 the thick films surrounding the 
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FIG. 4: (Colour online) The solvent density profiles 
A~ 3 pi(z,r) at the bulk state point X, around a pair of type 
A colloids of radius R = 4A _1 , separated a distance (a) 
h = 19.5A" 1 , (b) h = 16X' 1 and (c) h = 9. OA" 1 . The left 
hand plots are the density profiles for species 1 and the right 
hand plots for species 2. In (a) we see the two colloids are 
covered by thick wetting films rich in species 2, but they are 
far enough apart that the films do not connect. In (b) the 
colloids are sufficiently close that the films connect to form a 
bridge between the two colloids. In (c) the colloids are close 
to contact. The solvent mediated potentials corresponding to 
these profiles are denoted by □ in Fig. [6] 



pair of colloids join to form a bridge between the colloids 
composed of the wetting phase rich in particles of species 
2; see Fig. |4jb). For all separations h < hbn the bridge 
remains between the two colloids. In Fig. [6] we see that 
at h = h(, r , the separation at which the adsorbed films 
switch from the unbridged to the bridged state, the first 
derivative of W(h) changes discontinuously, i.e. there is a 
jump in the solvent mediated force / = — dW/d/iP^For 
h < hb r , W(h) decreases rapidly as h is decreased. We 
say more below about W(h) when there is bridging. 

In Fig. [5] we display the solvent density profiles around 
a pair of type B colloids immersed in the solvent phase 
X for the separations h = 9A _1 and h = 10A -1 . These 
'neutral' colloids are not covered by any thick wetting 
films. As the two B colloids are brought together, the 
solvent density profiles change continuously with h and 
we do not find any bridging. Comparing the density pro- 
files in Fig.[5^a) for h — 10A" 1 with those in Fig. [5jh) for 
h = 9A _1 , we see that as the pair of colloids are brought 
close together the solvent particles are expelled from the 
region directly between the two colloids. The SM po- 
tential W(h) corresponding to this case is displayed in 
Fig. [6] We find that W(h) is strongly attractive (near to 
contact, W(h = 9A" 1 ) = ~905k B T) and is fairly short 
ranged. When the two colloids are brought close to one 



£ 



£ 




P 2 (z,r) 
-i 1 r- 



OO 



_| I I 



- I 1 1 - 



CO 



_l I L_ 



-20 -10 10 20-20 -10 10 20 

Xz Xz 



FIG. 5: (Colour online) The solvent density profiles 
A -3 pi(z,r) at the bulk state point X, around a pair of type 
B colloids of radius R = 4A -3 , separated a distance (a) 
h = 10A -1 and (b) h = 9A _1 . The left hand figures are the 
density profiles for species 1 and the right hand figures for 
species 2. There is an adsorbed layer of both solvent species 
around the colloids but there are no thick wetting films. The 
solvent mediated potentials corresponding to these profiles are 
denoted by A in Fig. [6] 



another the expulsion of the solvent particles from the re- 
gion between them results in fewer solvent particles being 
directly at the surface of the colloids. This in turn means 
that the energy penalty for solvent particles being close 
to the colloids, resulting from the repulsive tails of the 
potentials <Pi C (r), is reduced. This lowering of the po- 
tential energy then leads to a lowering of the free energy 
of the system when the colloids are closer together and 
so W(h) is attractive. Below h ~ 9A _1 W(h) decreases 
extremely rapidly as the strongly repulsive parts of the 
colloid-solvent potentials start to overlap. 

In Fig. [6] we also display the solvent mediated poten- 
tial between a pair of type C colloids. Type C colloids, 
like type B colloids, are not covered by any thick wet- 
ting films. We find that the solvent mediated potential 
W(h) between a pair of type C colloids is almost indis- 
tinguishable from the same quantity for a pair of type B 
colloids for this state point X. The reason that the colloid 
B and colloid C solvent mediated potentials are so similar 
is that the two types of colloids differ only in the value 
chosen for the parameter 62 in Eq. (14), the magnitude 



of the interaction of the colloids with species 2 particles. 
At state point X, the density of the species 2 particles 
is small and so the potential energy contribution to the 
free energy from the interactions of the colloids with the 
species 2 solvent particles is also small. 

If we compare in Fig. [6] the solvent mediated poten- 
tial W(h) for the case when there are thick wetting films 
around the colloids (type A) and those when there are 
no thick films (type B or C colloids), it is tempting to ar- 
gue that by subtracting one potential from the other, one 
would be left with the contribution due solely to the thick 
wetting films present in the former case. We believe this 
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FIG. 6: The solvent mediated potential W(h), where h = 
zo — z\ is the distance between the centers of the colloids, for 
pairs of colloids of the same type, A, B, and C. The solvent is 
at state-point X in Fig.[T] For pairs of type B and C colloids, 
W(h) is a smooth continuous function of h and is almost iden- 
tical for both types of colloids. It is also strongly attractive: 
W(h = 9A _1 ) = -905fcgT. For type A colloids, W(h) is even 
more strongly attractive, W(h — 9A _1 ) = — 2442A; B T, and 
much longer ranged. We also observe a 'kink' (discontinuity 
in the first derivative) in W(h) at hbi — 18. 7A , where the 
thick films around the A colloids first connect to form a bridge. 
The symbols □ and A carry labels referring to Figs. [4] and 
[5] which display the solvent density profiles for the particular 
separation. 



argument to be cogent, at least at a qualitative level, on 
the basis of the following: Subtracting the value of the 
solvent mediated potential W(h c ) near contact at h = 
h c — 9A _1 for type C colloids from the same quantity for 
type A colloids gives us an estimate for the contribution 
to the solvent mediated potential due to the fluid bridge, 
i.e. W bl (h = h c ) ~ -2U0k B T + 900k B T = -1540fc B T. 
We now return to Eq. ( |10[ ) which is an expression ob- 
tained from macroscopic thermodynamic arguments for 
the contribution to W(h) arising from the presence of 
thick wetting films. In RefP^I the following approximate 
expression for the ratio S(h) of surface areas was obtained 
by assuming that the surface of the fluid bridge may be 
approximated by the surface generated by rotating the 
arc of a circle around the z-axis: 



S(h) 



wh(L 2 - (h/2 - w) 2 ) 1/2 (h/2-w 

arcsin 

4L(h/2 - w) 2 V L 



(L + h/2-w) 
2L 



w 



2L{h/2-wY 



(25) 



where L = R + I and 2w is the width along the z-axis 
of the bridge section - see RefP^ for further details, w 
is treated as a variational parameter, i.e. one selects the 
value of w that minimizes S(h). Using this expression 
for S(h), together with the value L = R + I = 7.5A -1 
(a rough estimate for the wetting film thickness based 
on inspecting the density profiles in Fig. |3ja)) and the 



value for the surface tension of the liquid-liquid inter- 
face 7 listed in Table [I] we obtain from Eq. (10 1 the 
value W(h c ) ~ — 1620/csT which is close to the value 
Wbr(h c ) — — 1540fcsT estimated above. Eqs. (10) and 



( 25 1 together also predict that bridging occurs when the 
colloids are at a distance h ~ 17. 2 A -1 , in fairly close 
agreement with the full DFT result of h^ r — 18.7A -1 . 
These observations lead us to conclude that the approxi- 
mation for W(h) in Eq. ( fTo| ), taken together with a good 
approximation for S(h), such as that in Eq. (25 1, provide 



a fairly reliable description of the solvent mediated poten- 
tial. Recall that Eq. (JTo|) was derived from macroscopic 
arguments, so it is perhaps surprising that the present 
microscopic DFT results for a system where the colloids 
are only one order of magnitude bigger in size than the 
solvent particles are well accounted for by this equation. 

The results presented in this subsection are for the col- 
loids immersed in the solvent at state-point X. However, 
similar results are obtained at state-point W which is 
slightly off bulk coexistence - see Table |T] Furthermore, 
the results obtained here for colloids of type A in the 
solvent at state point X are qualitatively the same as 
those for colloids of type C at state point Y, since in 
both cases the colloids are covered in thick films of the 
coexisting phase. Moreover the behavior of the systems 
with colloids of type A in the solvent at state point Y 
and colloids of type C in the solvent at state point X are 
also similar. The solvent around two colloids of type B 
behaves in much the same way in both bulk phases X 
and Y leading to similar solvent mediated potentials. 



D. Interaction between a single colloid and the 
fluid-fluid interface 



In Fig. [7] we display the density profiles for a single 
type A colloid at three different locations z§\ these lie on 
either side and within the fluid-fluid interface between 
coexisting phases X and Y. The plane of the fluid-fluid 
interface is set to be at z = with phase X being at z < 
and phase Y at z > as in Fig. |2ja). When the colloid is 
far from the interface at zq — > ±oo the colloid does not, 
of course, interact with the interface. Recall that col- 
loid A strongly prefers to be in solvent phase Y (see Sec. 
VB| above) and that when it is in phase X it is covered 
by a thick wetting film composed of the Y phase. This 
means that the global minimum of the grand potential 
for this system corresponds to zo — > oo, when the colloid 
is deep in the bulk phase Y. When the location of colloid 
A is shifted into the fluid-fluid interface the grand po- 
tential is minimized by maintaining a film rich in species 
2 particles around the colloid. This leads to bending of 
the fluid-fluid interface to accommodate the colloid - see 
Figs.[7ja) and^b). As z is further decreased, a point is 
reached where the grand potential cost from creating ad- 
ditional fluid-fluid interfacial area becomes equal to the 
grand potential for inserting the colloid into the bulk X 
phase. Beyond this point, the grand potential for the 
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FIG. 7: (Colour Online) The solvent density profiles 
X~ pi(z,r) around colloid A at locations: (a) zq — 0, (b) 
zo = — 20A" 1 and (c) zo — — 23A -1 in or near the fluid-fluid 
interface between co-existing phases X (z < 0) and Y (z > 0) 
in Fig. [l] The left hand density profiles are for solvent par- 
ticles of species 1 and the right hand profiles for species 2. 
The grand potential corresponding to these density profiles is 
denoted by □ in Fig. [9] 



bent interface configuration, such as that in Fig. |7jb) , is 
greater than that of having the colloid covered by a wet- 
ting film that is disconnected from the fluid-fluid inter- 
face; see Fig.[7jc). Decreasing z further does not change 
the grand potential. The SM potential (grand potential 
difference) AH,(zq) = Q,(zq) — ^( z o — > oo) is displayed 
in Fig. [9] We see that Afl(z a ) increases continuously as 
Zq is decreased until it reaches the point zq = — 22.1A -1 , 
where the film around the colloid detaches from the main 
fluid-fluid interface. This results in a discontinuity in the 
first derivative of AQ(z ) at this point. 

In Fig.|8]we display the density profiles for a single type 
B 'neutral' colloid at several different locations zq, i.e. 
either side and within the fluid-fluid interface between 
coexisting phases X and Y. The argument presented in 
Sec. |II C| applies in this situation. Due to the fact that 
the colloid does not have a strong preference for either 
solvent phase the global minimum of the free energy oc- 
curs when the fluid-fluid interface intersects the colloid; 
the density profiles for this situation are displayed in Fig. 
[8jc). If the colloid is moved to points on either side of 
the minimum at zq = 0.5A -1 , the fluid-fluid interface re- 
mains connected to the colloid and has an increased area. 
This increase in intcrfacial area results in an increase in 
the grand potential. The solvent density profiles for such 
configurations are displayed in Fig. |8jb) and[8jd). Fi- 
nally, as the colloid is moved even further from the in- 
terface a point is reached where the interface disconnects 
from the colloid. This is the case for the solvent density 
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FIG. 8: (Colour online) Same as Fig. [7]but now for colloid B 
at locations: (a) zq = 12A -1 , (b) zo = 9A _1 , (c) zo = 0.5A -1 
(the global minimum of the grand potential energy and the 
equilibrium position for the colloid), (d) zo = — 11A -1 and 
(e) zq = — 13 A -1 . The grand potential corresponding to these 
density profiles is denoted by A in Fig. [5] 



profiles displayed in Figs. |8ja) and |8je). In Fig. [9] we 
display Afi(zo) obtained from these density profiles. At 
the two values of Zq where the interface disconnects from 
the colloid, zq — — 11.7A" 1 and z$ = 9.5A -1 , we find a 
discontinuity in the first derivative of Af2(zo). Note that 
since we define AQ(zo) = fl(zo) — Q(zq — > 00), i.e. the 
grand potential when the colloid is at z minus that for 
placing the colloid deep into the bulk of the Y phase, 
we find that Af2 is zero for the density profiles in Fig. 
|8ja). For the density profiles in Fig. |8je), i.e. inserting 
the colloid into the coexisting X phase, Af2 = 403fcsT. 
That this value is greater than zero reflects the fact that 
the B colloids have slight preference for solvent phase Y 
over phase X. For the density profiles in Fig. |8jc), at 
the minimum of the grand potential, Af2 = — 312fcs7\ 
We may compare this result with the value one would 
obtain from the macroscopic thermodynamic approach 
discussed in Sec. |IIC| where we argued that in this sit- 
uation Aft is given by Eq. (Ill, i.e. roughly equal to 



the change in the area of the interface between the two 
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FIG. 9: The grand potential difference AQ(zo) = Cl(zo) — 
fi(oo) for a single type A or type B colloid in the vicinity of 
the fluid-fluid interface between co-existing phases X (z < 0) 
and Y (z > 0) in Fig. [I] zo is the location of the center of the 
colloid. We also display the solvent density profiles pi (z) and 
pi(z) for the unperturbed interface (Fig. |2|a)). The symbols 
□ and A carry labels referring to Figs.[7|or[8]which display the 
solvent density profiles around the colloid for that particular 
value of zq. We see that colloid A strongly prefers to be in 
phase Y (rich in species 2) and so the minimum of Afl(zo) 
for colloid A is at zo — > oo. The 'neutral' colloid B does 
not strongly prefer one phase over the other and there is a 
minimum in AQ(zo) at zq — 0.5A -1 , where the colloid sits 
within the interface; see Fig. |8|c). 
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FIG. 10: (Colour online) The solvent density profiles 
A pi(z, r) around a type A colloid located near a planar wall 
at z — 0, with potentials given by Eq. (15 1. In (a) the centre 
of the colloid is at zo — 21A -1 , and the colloid is sufficiently 
far from the wall that the adsorbed films do not interact. In 
(b) zo — 18A -1 and the wetting films covering the wall and 
the colloid connect. In (c) Zo = 5A _1 , and the colloid is close 
to contact with the wall. The bulk solvent is at state point 



W in Fig. [T] and is rich in species 1. The left hand profiles are 
for species 1 and the right hand profiles for species 2. The 
grand potential for these density profiles is denoted by □ in 
Fig. fT3l 



coexisting solvent phases when the colloid sits at the in- 
terface, —ttR 2 , multiplied by the fluid-fluid surface ten- 
sion 7. Using the value for 7 given in Table |T] we find 
Ail = -irR 2 j = -270k B T. One could also argue that 
because of the large values chosen for the parameters bi 
in Eq. (14), the effective radius R* of the colloid is some- 



what larger than R. From inspecting the density profiles 
in Fig. [3] we estimated R* ~ 4.5 A -1 (recall R = 4A _1 ), 
giving an estimate for Af2 ~ ~ttR* 2 ^ — —350ksT, which 
is in reasonable agreement with the value we obtain from 
the microscopic DFT theory. What this shows is that one 



may use the crude approximation in Eq. (Ill to estimate 



roughly Aft for 'neutral' colloids, even when the colloids 
are only one order of magnitude larger in size than the 
solvent particles. 



Colloid C strongly prefers solvent phase X so when we 
insert colloid C into the fluid-fluid interface it behaves in 
the opposite way to colloid A, i.e. in the phase rich in 
species 2 (state-point Y) colloid C is covered by a thick 
layer rich in species 1. This means that the shape of the 
curve AH(zq) for colloid C (not displayed) is similar to 
the curve for colloid A in Fig. [9] but with the replacement 

Zo -> -Zq. 



E. Interaction between a single colloid and a 
wetting film adsorbed at a planar wall 

In Fig.[lO]we display the solvent density profiles around 
a single type A colloid positioned at three different dis- 
tances zq from a planar wall, located at z = 0, with 
potentials given by Eq. ( [T5| . As the bulk solvent is at 
state point W in Fig. [T] which is near to phase coex- 
istence, the wall is covered by a thick wetting film of 
the coexisting phase rich in species 2. In the bulk fluid, 
away from the wall, colloid A is also covered by a thick 
wetting film. When the colloid is far from the wall, 
the films do not interact - see for example the density 
profiles in Fig. [To|a). As the colloid is brought closer 
to the wall the solvent density profiles change abruptly 
from a configuration such as that displayed in Fig. 10 'a) 



to a configuration where there is a bridge of the wet- 
ting phase connecting the colloid to the wall; see Fig. 
10b). The position of the colloid when the bridge forms 
is z = z t = 19.8A -1 . As zq is decreased further, the col- 
loid remains within the wetting film all the way to contact 
with the wall at zq = 4A _1 . In Fig. 13 we display the 
solvent mediated potential (grand potential difference) 
AH(zq) = il(z n ) — il(z — > 00) calculated from these den- 
sity profiles. We see that for zq < zt, AH(zq) is negative 
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FIG. 11: (Colour online) Same as Fig. 10 but for a type B 
colloid. In (a) zq = 14A" 1 and the colloid is sufficiently far 
from the wetting film at the wall that these do not interact. 
In (b) zq — 10 A -1 and the film at the wall bends so that the 
interface intersects the colloid. In (c) zq — 5.2A -1 and the 
colloid is nearly in contact with the wall. The grand potential 
corresponding to these density profiles is denoted by A in Fig. 
[131 



indicating that the effective interaction between colloid 
A and the wall is strongly attractive. When the colloid 
is close to contact with the wall Afl(zo) — — 7050fcsT. 
The primary origin of this attraction lies in the fact that 
when z < z t the area of the interface between the bulk 
fluid and the wetting film covering the wall and the col- 
loid is less than when the colloid is in the bulk away from 
the wall at z$ > z t . The abrupt change in the density 
profiles at z a = z t manifests as a discontinuity in the first 
derivative of ACl(zo), i.e. a jump in the SM force between 
the colloid and the wall at zq = z t . 

In Fig.[TT]we display the solvent density profiles around 
a single type B colloid positioned at three different dis- 
tances zq from the wall. Colloid B is not covered by a 
wetting film in the bulk and as it is brought closer to 
the wall there is an abrupt change in the solvent den- 
sity profiles at zq = z t = 13.3A -1 , from a configuration 
where the wetting film covering the wall is almost com- 
pletely unperturbed, to a configuration where the wetting 
film extends to meet the colloid, so that the interface be- 
tween the bulk fluid phase and the wetting film is an- 
chored to the colloid. As zq is decreased further this 
interface remains connected to the colloid; see Fig. [Tl|b) 
and (c). In Fig. 13 we display the grand potential differ- 
ence Afi(zo) = Q(zo) — Q( z o —> °°) calculated from these 
density profiles. We see that at zq — z t there is a discon- 



tinuity in the first derivative of AQ(zq), and that Af2(zo) 
is negative for zq < z tl indicating that the effective in- 
teraction between colloid B and the wall is attractive. 
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FIG. 12: (Colour online) Same as Fig. 

colloid. In (a) zq — 8A~ X . Since the colloid strongly prefers 
to be in the bulk solvent phase the wetting film at the wall 
thins so that the colloid can maintain a layer of the bulk phase 
around itself as it approaches the wall. In (b) zq — 4.5A -1 
and the colloid is almost in contact with the wall. The grand 
potential corresponding to these density profiles is denoted by 
o in Fig. [13] 
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FIG. 13: The grand potential difference AQ,(zq) = £l(zo) — 
fi(oo) (solvent mediated potential) for a single type A, B, or 
C colloid, zq is the distance of the centre of the colloid from 
a planar wall at z = with potentials given by Eq. 1 15 1. We 
also display the solvent density profiles pi{z) near the wall 
when no colloid is present (Fig. [2jb)) ; the wall is covered by 
a thick wetting film rich in species 2. Symbols □, A and O 
carry labels referring to Figs.[l0]-ri2]which display the solvent 
density profiles around the colloid for these particular values 
of zq. The bulk solvent is at state point W in Fig. [T] and the 
vertical line is at zq = 4A _1 , the colloid radius. 



Close to contact AH,(zq) ~ — 1500fcsT, so the effective 
interaction between colloid B and the wall is much less 
strongly attractive than that between colloid A and the 
same wall. 

In Fig.[l2]we display the solvent density profiles around 
a single type C colloid positioned at two different dis- 
tances zq from the wall. Colloid C strongly prefers to be 
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in the bulk solvent phase, which is poor species 2, rather 
than in the phase wetting the wall, which is rich in species 
1. As colloid C is brought close to the wall it moves into 
the vicinity of the wetting film covering the wall and be- 
cause colloid C prefers the bulk phase, the wetting film 
at the wall is thinned (and eventually expelled) in the 
region between the colloid and the wall. The density 
profiles change continuously as zq is varied. In Fig. [13] 
we display the grand potential difference Af2(zo). As 
zq is decreased AH(zq) increases smoothly as the grand 
potential for having the bulk solvent phase close to the 
wall is greater than for having the wall covered by a uni- 
form thick wetting film of the coexisting phase. The fact 
that near the wall ASl(zo) is positive indicates that the 
effective interaction between colloid C and the wall is 
repulsive, in marked contrast to colloids A and B. 



VI. DISCUSSION AND CONCLUDING 
REMARKS 

We have investigated the behavior of model colloids in- 
serted into a binary mixture of softcore solvent particles 
that exhibits fluid-fluid dcmixing. We considered solvent 
statepoints at, or close to, fluid-fluid coexistence where 
the solvent correlation functions are reasonably well de- 
scribed by a simple RPA DFT. By varying the colloid- 
solvent interaction parameters we were able to investigate 
the colloid's affinity for insertion into one of the bulk co- 
existing phases. For type A and C colloids we showed 
that when these are in isolation in the unfavored bulk 
solvent phase they adsorb a thick 'wetting' film of the 
favored coexisting phase. Type B colloid has no strong 
preference for either bulk phase. 

Using a 'brute-force' DFT approach that models the 
colloids by an external potential acting on the solvent 
particles we calculated the density profiles of the solvent 
for the following scenarios: (i) a pair of colloids of the 
same type immersed in the bulk, (ii) a colloid in a pla- 
nar fluid-fluid interface, and (iii) a colloid near a planar 
wall that adsorbs a thick wetting film. By calculating 
the grand potential of the solvent we were able to de- 
termine the SM interaction in each case. We found that 
the presence (or absence) of a thick adsorbed film around 
the colloids in isolation determines the behaviour of the 
system in the different scenarios. 

When two colloids, both with a thick adsorbed film, 
are brought together then the films can connect together 
abruptly to form a fluid bridge. A similar phenomenon 
occurs when a colloid is brought close to a wall. The 
formation of a fluid bridge gives rise to a strong, long- 
ranged SM attraction. If there are no thick wetting films 
then the SM potential is much less attractive and much 
shorter ranged. Furthermore, when the wall and colloid 
have very different solvent affinities then the SM interac- 
tion can be repulsive. 

For a type B colloid inserted in the fluid-fluid interface 
the effective colloid- interface interaction has an attractive 



well implying that such 'neutral' colloids would become 
trapped in the interface. On the other hand, colloids 
A or C that strongly prefer to be in one of the bulk 
phases distort the interface, and ultimately experience a 
repulsive interaction with the interface. 

In our mean-field DFT treatment we find an abrupt, 
first-order like transition for many of the scenarios stud- 
ied, i.e. the SM force is discontinuous when a fluid 
bridge forms. However, since a finite number of par- 
ticles arc involved in forming the bridging film and in 
the onset of interface distortion, these cannot be true 
phase transitions. In reality fluctuation (finite size) ef- 
fects must result in a rounding of the discontinuity in the 
SM force P2 The authors of Refill suggested that a crude 
estimate for the width of this rounding, Szq ^, can be 
obtained by considering that fluctuations should only be 
relevant when | An Brid s cd (z ) - Aft Unbrid s cd (z )| < k B T 
where Aft (Un)Bridgcd (z ) is the grand potential on the 
(un)bridged side of the transition. Using this crude cri- 
terion we find that the transition is smeared over a length 
^o.br/zo.br ~ 10~ 3 . This is similar to the value estimated 
by Archer et alW^ in their DFT studies of bridging be- 
tween two big solute particles in a Gaussian core model 
of a binary solvent. 

We showed that for two colloids in bulk the strength 
and range of the SM interaction is reasonably well de- 
scribed by a simple (capillarity) model for the shape of 
the liquid bridge formed between the two colloids, namely 
Eqs. ([9| and (251. The only inputs to this theory are the 
surface tension of the planar interface between the two 
solvent phases and the thickness of the adsorbed 'wetting' 
film. These quantities can be obtained from simpler DFT 
calculations since they depend only on profiles that vary 
in one dimension. For the neutral colloid immersed in 
the planar fluid-fluid interface we calculated the depth 
of the attractive well in the effective colloid-interface in- 
teraction and found this to be reasonably close to the 
value predicted from a very simple model that considered 
only the surface of the colloid in an unbending interface, 
Eq. (111. It is interesting that these macroscopic (cap- 
illarity) approximations appear to have some validity in 
the case where the colloids are only one order of mag- 
nitude larger than the solvent particles. Furthermore, 
more systematic comparisons of the results of full DFT 
calculations with those based on macroscopic approaches 
would be valuable in ascertaining the limitations of the 
latter for nanoparticles. 

We have deliberately limited the scope of this study 
and have not conducted an extensive investigation of the 
entire phase diagram, nor of the full set of parameters 
characterizing the solvent-colloid and solvent-wall inter- 
actions. It is likely that the presence of a pre- wetting 
transition would lead to additional features in the SM 
potentials. Furthermore, we have only investigated prob- 
lems where the density profiles exhibit cylindrical sym- 
metry. The interactions between three hard-sphere col- 
loids in a bulk hard-sphere solvent have been investigated 
both with brute force DFT 51 , and the particle insertion 
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methods, and an obvious extension of the present work 
would be to consider the interactions between three col- 
loids in the bulk, or two colloids in the presence of a fluid 
interface, or near a planar wall. This would necessarily 
require a full three dimensional DFT computation, which 
in turns brings its own complications. The density pro- 
files can no longer be easily calculated with sufficient pre- 
cision on a desktop computer but instead one might use a 
number of computers in parallel.^ Furthermore, partic- 
ular care must be taken in establishing that the density 
profiles are the true equilibrium profiles and that any 
abrupt jumps correspond to the equilibrium phase tran- 
sitions, i.e. the grand potential must be calculated very 
precisely; this poses considerable numerical challenges. 

Our approach can be compared to investigations us- 
ing both molecular dynamics and Monte Carlo simula- 
tions. Many of the studies mentioned in the Introduction 
used Lennard- Jones or similar potentials for the solvent- 
solvent and solvent-colloid interactions, and further work 
based on the DFT approach should attempt to treat such 
models. On the other hand, a particular advantage of the 
present simple model is that the ratio between the col- 
loid and solvent particle sizes is sufficiently large that we 
are able to obtain a thick film adsorbed around modest 
sized (soft core) colloids. Existing simulation studies us- 
ing Lennard- Jones potentials have used a maximum col- 
loid diameter of 10 times the solvent particle diameterPS 
This is not sufficient to adsorb a thick wetting film for 
hard core colloids. 



So far in this paper we have not indicated the form of 
the bare interaction between the large colloids, nor the 
interaction between the colloids and the wall, since these 
have no bearing on the SM potentials. We suggest that 
the bare colloid-colloid potential has a hard-core of diam- 
eter 8 A -1 and that the potential decreases rapidly out- 
side this hard-core. Then the effect of the bridging film is 
to still induce a strong, long-ranged attraction between 
the colloids for suitable solvent statepoints. Similarly, 
the bare wall-colloid potential should diverge rapidly at 
z = 4A~ X , and for suitable solvent state-points the long- 
ranged attraction should be retained. Finally, although 
our model does not correspond directly to an experimen- 
tal situation it has been suggested on general grounds 
that light scattering experiments could be used to in- 
vestigate the value of the second virial coefficient, B2 of 
colloids in bulk solvents.^ Rapid changes of B2 upon 
changing the solvent state-point might signal the aggre- 
gation of colloids driven by the formation of fluid bridges. 
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Note that in Eq. (|2| we have implicitly assumed that the 
colloid is hard (i.e. the potentials 4>t c (r) that it exerts on 
the solvent particles have the property 4>T{r < R) = oo) 
and therefore choosing the Gibbs dividing surface at r = R 
is the natural choice. One can, of course, define the dividing 
surface at any value of r without changing Q,. The choice of 
location for the dividing surface does, however, affect the 
values for the available volume and 7 pa used in Eq. j2J|. 
In order to to evaluate W(h) it is necessary to calculate 
the series of solvent density profiles for various separations 
h and it is advantageous to use the density profile solu- 
tion obtained from a neighboring value of h as an initial 
estimate. If we start with the colloids located far apart 
and decrease h then we calculate a branch of W(h) cor- 
responding to the non-bridged configuration. If instead we 
start with the colloids in close proximity and increase h 
then we calculate another branch of W(h) corresponding 
to a bridging film between the colloids. Where these two 



branches cross, i.e. where the grand potentials are equal at 
h — hbi, then the two solutions coexist, and there is a dis- 
continuity in the SM force / = —AW /Ah. The non-bridged 
branch of W(h) is metastable for h < hhr and the bridged 
branch is metastable for h > h^ r ; more details are given in 
RefP. 



